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Abstract. Recent experimental data on dijet cross sections in diffractive photoproduction at HERA collider 
are analysed with an emphasis on QCD factorisation breaking effects. The possible sources of the different 
conclusions of Ifl and ZEUS collaborations are studied. 
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The first observation of deep inelastic scattering (DIS) 
events at HERA containing a large gap in the pseudora- 
pidity distributions of final hadrons [Tl[2] has generated a 
considerable renewed interest in understanding of diffrac- 
tion. In the ep diffractive interactions measured at the 
HERA experiments the proton can stay intact or dissoci- 
ate into a low mass state (Y) , while the photon dissociates 
into a hadronic state X, 7*p Xp'{Y). The systems are 
separated by a non-exponentially suppressed large rapid- 
ity gap being present due to vacuum quantum numbers 
of the diffractive exchange [3]. The diffractive exchange, 
called pomeron (IP), carries away a fraction xp of the ini- 
tial proton longitudinal momentum and a four-momentum 
transfer t = {P — PyY associated with the outgoing (lead- 
ing) system yQ. 

One of the most interesting questions which are dis- 
cussed in the studies of diffractive processes is whether 
they can be considered being factorisable into diffractive 
parton distribution functions (DPDFs) of the proton and 
perturbatively calculable partonic cross sections. Such a 
concept, called hard QCD (or collinear) factorisation, was 
proved to be valid in the regime of diffractive deep inelastic 
scattering (DDIS) (4j . On the other hand the factorisation 
concept was found not to be valid for hard processes in 
diffractive hadron-hadron scattering as measured for ex- 
ample in pp interactions at Tevatron . The predictions 
of cross sections of diffractive dijet production in pp based 
on DPDFs provided by the HI collaboration, overestimate 
the measured cross section almost by one order of mag- 
nitude. This observation was explained later on theoreti- 
cally [6j assuming that factorisation breaking results from 
absorptive effects caused by multiple rescattering effects, 
see e.g. [71IH| for earlier discussion. The rapidity gap can 



where P is the four-momentum of incoming proton 



then be populated by secondary particles which spoil the 
experimental signature of the diffractive event. 

At HERA, diffraction has been studied in a wide range 
of exchanged photon virtualities, Q^, ranging from photo- 
production regime, ~ 0, to the deep inelastic scatter- 
ing with ^ 0. Various hard diffractive processes were 
studied at HERA such as inclusive ones, production of 
jets, charm and vector mesons. 

The factorisation in diffractive DIS dijet production 
was experimentally tested by the HI and ZEUS collabo- 
rations pillO). Far from clear, however, is the situation in 
the photoproduction regime of the ep diffractive scatter- 
ing. In photoproduction, in the leading order approxima- 
tion, the small photon virtuality allows for partonic fiuctu- 
ations that live long enough. The photon may not couple 
directly to the quarks in the proton, but only a part of its 
four-momentum participates in the hard interaction. Such 
interactions are called resolved. The photon can still cou- 
ple directly (with its whole four-momentum) to the quarks 
and these interactions are called direct. The resolved pho- 
ton interactions resemble the hadron-hadron ones since 
two particles with structure scatter on each other. The 
variable x~f, which is defined as a four-momentum frac- 
tion of the photon taking part in the hard interaction, is 
used to distinguish between the two regimes in photopro- 
duction. Obviously, following relations hold: a;^ = 1 and 
x-y < \ for the direct and resolved photon interactions, 
respectively. Effects of fragmentation and a finite experi- 
mental resolution impose a smearing on the value of x-y . In 
experiments, often a value of x~f around 0.75 is considered 
to be a discriminator below (above) which the events are 
regarded as being due to resolved (direct) photon inter- 
action, with reasonably low contamination of one in each 
other. 

Because of the presence of these two photon interaction 
regimes, studying diffractive processes in photoproduction 
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Table 1. Kinematic cuts used in HI ^ and ZEUS [18] 

analyses 



HI 


ZEUS 


Q2 < 0.01 GeV^ 


< 1 GeV^ 


0.3 < 2/ < 0.65 


0.2 <y < 0.85 


> 5 GeV 


> 7.5 GeV 


> 4 GeV 


> 6.5 GeV 


-1 < 77J''"(2) ^ 2 


-1.5 < ryj""*^' < 1.5 


DifTractive cuts 


xp < 0.03 


Xp < 0.025 


zp < 0.8 




\t\ < IGeV^ 


|t| < IGeV^ 


My < 1.6 GeV 


My = mp 



is a useful tool for tests of the validity of the hard factori- 
sation in diffraction. The partonic structure of the photon 
is described in terms of quark and gluon densities that 
obey DGLAP evolution equations [TT]. There was an ear- 
lier prediction of Kaidalov et al. [H] that the factorisation 
breaking of the resolved part should induce a suppression 
of the NLO QCD expectation by about a factor of 0.34. 
This idea was widely discussed and applied to published 
data in the studies of Klasen and Kramer [TSIIT^ . However 
recently [TS] theoretical expectations were revised stress- 
ing the fact that due to the inhomogeneous term in the 
DGLAP evolution there is also point-like part of the pho- 
ton structure function [16] . The hadron-like part of the 
photon structure (suppressed by 0.34) occurs only at low- 
est values of four-momentum fractions Xj ^ 0.1 which 
are experimentally hardly accessible. The dominant mea- 
sured part of the resolved processes is therefore induced 
by the point-like component of the photon structure func- 
tion with a significantly weaker suppression as compared 
to the 0.34, see [15J for details. 

2 Recent results from HERA 

The HI and ZEUS collaborations analysed the diffractive 
dijet photoproduction data in [TTlfTSlfT^ . The advantage 
of using the diffractive photoproduction of dijets is that 
two jets in the final state enable us to reconstruct the 
variable. The HI collaboration observed a global suppres- 
sion of the dijet cross sections with respect to NLO QCD 
calculations by a factor of about 0.5 [T7]. On the contrary 
the ZEUS data published in [18,19^ were compatible with 
no suppression. A new study of HI collaboration [20j with 
a three times larger data sample and a similar kinematic 
region as the previous study [17\ fully confirms the previ- 
ous HI observation. In contradiction with expectations of 
only the resolved processes being suppressed [T^[T5] nei- 
ther experiment did observe a difference in suppression 
for the resolved enriched part {x^ < 0.75) and the direct 
enriched part of the cross section {x~^ > 0.75). 

In HI and ZEUS analyses the diffractive events were 
selected with a large rapidity gap method and jets were 
identified using the inclusive /ct cluster algorithm |21) in 
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Fig. 1. The diffractive quark singlet and gluon densities for 
the squared factorisation scale fi^ — 25 GeV^ in the region 
My < 1.6 GeV and [i| < 1 GeV^ for two values of xp, xp — 
0.01 (left) and xp = 0.02 (right). Uncertainties of the HI fit B 
are depicted as a band. 

the laboratory frame. The phase space of both analyses 
was different (see Table 1), the main difference being at 
somewhat larger transverse momenta of the jets, Ei^^ ^ of 
ZEUS measurement. HI collaboration collected data with 
a tagged electron that allowed to restrict the to very 
low values (Q^ < 0.01 GeV^), the ZEUS analysis was 
done with an untagged electron sample of events {Q^ < 
1 GeV^). In HI analysis the additional cut zp < 0.8 was 
applied since the DPDF sets are not valid at the largest 
values of zp 

The NLO QCD calculations were performed by means 
of using the program of Frixione et al. [5^ (HI) and Klasen 
and Kramer [23 (ZEUS). Both calculations provide a con- 
sistent results as it was demonstrated in j^ . 

For the NLO QCD calculations, the HI collaboration 
used three DPDF sets, HI 2006 DPDF fit B, measured 
in the analysis of the inclusive DIS data '5S^, the inclu- 
sive and dijets combined fit HI 2007 DPDF fit Jets [9] 
and inclusive and dijets combined ZEUS fit SJ [T^. HI 
and ZEUS collaborations compared their results with the 
calculations using the same DPDF fits. 

In Fig. 1. the comparison of the three DPDF sets is 
shown as a function of the partonic longitudinal four- 
momentum fraction with respect to the pomeron, z, for 
two values of xp. As expected a good agreement of all 
three fits is seen for the quark singlet density. The dif- 
ferences in the fits are much larger for gluon densities 
which are more important for dijet measurements. Note 
that the values of ZEUS fit SJ in the Fig. 1 are multi- 
plied by a factor of 1.2. This is connected with the differ- 
ence in methods which were used to measure diffractive 
inclusive cross section by HI and ZEUS collaborations. 
The selection of diffractive events relying on the rapid- 
ity gap method yields a sample which is dominated by 
elastically scattered protons, but which also contains an 

^ Zp is the longitudinal four-momentum fraction of the par- 
ton entering the hard subprocess with respect to pomeron. 
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Fig. 2. Differential cross sections for diffractive dijet photoproduction as published by ZEUS [18j compared with NLO QCD 
calculations corrected for hadronisation using DPDF ZEUS fit SJ (white line), and two Ifl DPDFs divided by 1.2 - HI fit B 
(full line) and HI fit Jets (dashed line) . For NLO QCD calculations ZEUS fit SJ, the theoretical uncertainties connected with 
the change of scale are shown as a dark band. 



admixture of events in which the proton dissociates to low 
mass {My) state. In the case of HI collaboration the mea- 
surements are corrected to the region My < 1.6 GeV and 
|t| < 1 GeV^ while ZEUS collaboration corrects the mea- 
surements to My — nip, (for more details see e.g. [25]). 

3 Crosscheck of HI and ZEUS results 

In the following the crosscheck and discussion of the possi- 
ble sources for different HI and ZEUS conclusions as well 
as the sensitivity to an alternative photon parton distri- 
bution function (7-PDF) choice will be studied. 

Hadronisation corrections - To be able to compare the 
NLO QCD cross sections calculated at the level of par- 
tons with the data, the correction for hadronisation effects 
must be done. A common practice is to use Monte Carlo 
leading order models with initial and final state parton 
radiation in order to mimic some features of the next-to- 
leading order parton dynamics. In |15] the justification of 
this procedure is discussed. It is noticed that the hadro- 
nisation corrections based on the LO Monte Carlo should 
be different from the effects associated with NLO calcula- 
tions. To minimize this effect both collaborations used the 



procedure of reweighting of MC parton level distributions 
to NLO calculated parton level distributions. Such a pro- 
cedure ensures that the bin-to-bin migrations due to the 
transitions from partons to hadrons will be described more 
accurately. The hadronisation corrections calculated and 
applied to NLO predictions were however different in HI 
and ZEUS analyses. The hadronisation corrections, Shadr, 
(defined in each bin i as (1 -I- Shadr) i = Ci — crf°'^''/crf"'^*, 
where Ci is the correction factor to be applied to the NLO 
prediction) obtained in ,17 and independently also in |20] 
are on average at the level of —15%. The hadronisation 
corrections calculated in [T^ and applied in [TSlfT^ are 
also at the level of 15% but positive. Nevertheless, it is 
evident that the difference in hadronisation corrections 
cannot explain the contradiction in conlusions of the two 
analyses. Identical hadronisation corrections would even- 
tually make the difference in suppression even larger. 



An attempt was made to recalculate the hadronisa- 
tion corrections in the kinematics of the analyses of both 
HI and ZEUS experiments. The recalculated hadronisa- 
tion corrections for HI kinematics agree within few per- 
cent with the published |5D] ones while the recalculated 
hadronisation correction factors for ZEUS kinematics are 
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Fig. 3. The differential cross sections for diffractive dijet photoproduction as measured by HI :20, compared with NLO 
calculations using HI fit B DPDF and photon structure functions GRV and AFG. The dark bands around the theoretical GRV 
PDF predictions (marked by white line) show the uncertainty of NLO calculations connected with the scale variation by factors 
0.5 and 2.0. 



lower by about 15% than those from [T^. In Fig. 2, the 
measured ZEUS cross sections for variables x^, i?;^*^, xp, 
Zp, invariant mass of diffractive hadronic system Mx and 
inelasticity y from [18] are compared with the NLO pre- 
dictions corrected by means of using of the recalculated 
hadronisation corrections. The experimental differential 
cross sections are clearly in agreement with NLO calcula- 
tions using all DPDFs except for da/dzp where the shapes 
of experimental and predicted distributions disagree. As 
expected the conclusions of both analyses as concerns the 
non-observation or observation of the factorisation break- 
ing are not sensitive to hadronisation corrections used. 



Alternative photon distribution function - In the analy- 
ses mentioned above the photon GRV structure function 
P7] was used in the NLO calculations. It was noticed jT5] 
that the point-like part of the resolved contribution in 
GRV 7-PDF may be overestimated by about 25% in com- 
parison with the more recent AFG 7-PDF [55]. In Fig. 3 
the HI experimental differential cross sections |20] are 

^ Note that hadronisation corrections in ZEUS case were 
evaluated in [18] using the older DPDF than DPDF ZEUS 
fit S.] used here. 



shown for variables: x^, Ei^ , log^^Q xp, zp, mean pseudo- 
rapidity of two jets, {rj^'^^^), absolute value of the difference 
in pseudorapidities of both jets, |Z\?7j°**'|, total hadronic en- 
ergy W, the invariant mass of jets M12 and Mx- These 
experimental distributions are compared with NLO QCD 
predictions based on GRV and alternatively the AFG 7- 
PDF corrected for the effects of hadronisation. Indeed, 
it can be seen that the NLO QCD calculation based on 
AFG 7-PDF predicts smaller cross section than the one 
provided by GRV 7-PDF. In any case the difference in 
the value of suppression is not significant and is covered 
by the experimental uncertainties and by the variation of 
the renormalization and factorisation scale (yUr j/2, 2/Zr,/). 
We can summarize that the conclusions of both HERA ex- 
periments are not sensitive to the 7-PDF used. 



IVIatching the different phase space - The HI and ZEUS 
dijet data cannot be compared directly, since they have 
different kinematic domains. It was suspected that con- 
tradictory conclusion on the values of suppression is due 
to the different cuts on the of the jets in both analyses. 
This was not fully confirmed, a tendency of weakening of 
the suppression with Et increasing is observed but within 
the uncertainties it is still present [5^150] . Here an at- 
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Fig. 4. Differential cross sections for diffractive dijet photoproduction as published by ZEUS [18] (points and full line) compared 
to extrapolated Ell differential cross section (squares and dotted line). 



tempt is made to recalculate (using Monte Carlo model 
RAPGAP [26]) the results from HI kinematic domain to 
the one of the ZEUS analysis, see Table 1. 

Conversion from the HI phase space to that of ZEUS 
could be done in two steps. 

In the first step the differential cross sections measured 
by HI are multiplied in each bin by factors, obtained by 
MC RAPGAP, accounting on the extension from < 
0.01 GeV^ to Q2 < iGeV^. This leads to an increase of 
the differential cross sections by about a factor of 1.35 with 
no substantial change of their shapes. We expect that the 
uncertainty of the extrapolation in can be neglected 
since it is given by quantum electrodynamics only - in 
these analyses is used as a hard scale. Next, the data 
are divided by a factor of 1.2 to correct the cross section 
to the elastic case of My = mp of ZEUS0 

In the second step a matrix is constructed, by means of 
MC RAPGAP, which accounts on transition from the HI 
phase space to that one of ZEUS, 'M.h^z- The matrix is 
determined separately for each pair of corresponding dif- 
ferential cross sections from both analyses (x-y, i?^*^,. . . ). 
The matrix contains the probabilities that a particular 
event belonging to the i-th bin will end up in the j-th 
one of the given cross section as one makes the transition 
from the HI to the ZEUS analysis phase space. Although 
the shapes of differential cross sections measured by HI 
are rather well described by MC [2D], to achieve a higher 
accuracy of the matrix determination, the MC spectra are 

The ratio 1.20 ± 0.11 is taken from [29] . 



reweighted to match the data better. In addition the con- 
tribution of events generated outside the HI kinematics 
but fulfilling ZEUS cuts needs to be taken into account. 
The normalization of this contribution is provided by the 
previously discussed reweighting of MC. The including of 
these events depends on the validity of used model in the 
region uncovered by HI data. For each variable the result 
of the extrapolation procedure is given by a histogram az 
determined as 

az = Mh^z cjh + (jT (1) 

where an is the vector of values of the measured cross 
sections modified in step one and cr|'*'' is the contribution 
from outside of the HI phase space. 

Results of this transformation of HI data are presented 
in Fig. 4 together with ZEUS published data. It can be 
seen that after the transformation to the identical phase 
space the HI differential cross sections are lower than the 
ZEUS results by factor about 0.6. On the other hand the 
shapes of HI and ZEUS differential cross sections are very 
similar, except for zp. The uncertainties of the extrap- 
olated HI data were determined by propagation of the 
statistical and uncorrelated systematic errors from |20) . 
The relative uncertainty of a'^^^ was assumed to be the 
same as the one of the total cross section from |5D1. In 
addition the model dependence of the extrapolation was 
studied by using two different DPDFs - HI fit B and HI 
fit Jets. This source of uncertainty was found to be rela- 
tively small (with the exception oi zp distribution) and is 
included into errors bars shown in Fig. 4. 
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It is worthwhile to mention that two independent cross- 
checks of the transformation method were done. The first 
one is based on the comparison of HI resuhs from pp] 
extrapolated to the kinematic region of HI preliminary 
high E^^ data [24] . The extrapolated results agree within 
errors with cross sections from j24| in the shape and nor- 
malisation very well. The another crosscheck used ZEUS 
published data [T5] transformed to the HI preliminary 
high E^r^^ data [H] yielding the same relative difference 
as shown in Fig. 4. Note that in this case the HI cuts 
were the subset of ZEUS cuts and therefore the constant 
term of the equation (1) was not needed. 



4 Summary 

The analysis of published HI and ZEUS results is done 
with the emphasis to understand better the possible sources 
of the discrepancies between the conclusions of both col- 
laborations. It is shown that the results are not signif- 
icantly sensitive to the photon structure function used. 
Although, here the different hadronisation corrections are 
obtained than in |18l . it has no impact on the interpre- 
tation of ZEUS differential cross sections. The conversion 
of HI results to ZEUS phase space is done. The shapes 
of differential cross sections measured by both collabora- 
tions are in agreement (except for zp variable), however, 
the HI results are on average lower by about 40% than 
ZEUS ones. Within the limitations of the transformation 
method (based on Monte Carlo) there is a suggestion that 
the observed discrepancy between HI and ZEUS results 
concerning factorisation breaking is not caused by differ- 
ent phase space of both analyses. 

The puzzle of factorisation breaking in diffractive di- 
jet photoproduction could be resolved therefore by new 
experimental analyses. Among them the most promising 
one could be the identification of diffractive events based 
on leading proton detection. This method has both an 
advantage of providing a data sample free of proton disso- 
ciation and of reducing the uncertainties of the diffractive 
selection if compared with large rapidity gap method. 
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